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1. Introduction 



For many years there exists a permanent non passing away interest to the nu- 
cleon spin problem: how the miclcon spin is distributed among its constituents - 
quarks and gluons. This interest arose essentiaUy after appearance of the results 
of famous the EMC experiment ^ on deep inelastic scattering (DIS) of longitudinal 
polarized muons on longitudinal polarized protons. EMC, using also the earlier 
data of SLAC ^, came to the surprising and contraintuitive conclusion that quarks 
are carrying a small part of the proton spin projection in the polarized proton. 
This statement results in excitement of one part of particle physics community and 
to disappointment of the other. Many ideas, how to overcome this problem were 
suggested, sometimes new and interesting, sometimes, however, misleading. As a 
result of these investigations we understand now much more not only about pro- 
ton spin structure, but also about some others connected with its issues, and this 
understanding will be with us forever. 

In my talk I will try to lay stress on this new understanding, which arised in 
the last years, taking in mind the aphorism: "The result of calculation is not the 
number, but the understanding ". 

2. General Relations * 

The polarized deep inelastic e(yu)-nucleon scattering cross section is determined 
by imaginary part of the forward virtual photon-nucleon scattering amplitude 
T^vip-il)-! antisymmetric in i/. The general form of ImT^j^{p,q), following from 
gauge, P and T invariance is 

27r 



^ m 



(1) 



where p and q are the nucleon and virtual photon momenta, s and m - are the 
nucleon spin and mass. The structure functions Gi and G2 depend on two invariants, 
(f — —Q^ < and u = pq. When and u are large enough, but their ratio 
X — Q'^/2u is fixed, < x < 1 (so called scaling limit), the scaling relations take 

{u/m')G^{u,q^) = g,{x,Q') 

{v/m'^fG.M) = g2{x,Q^) 
In the parton model gi and g2 are functions of only scaling variable a;, in QCD 
a smooth (logarithmic) dependence of (71,(72 on appears. In the parton model 
gi{x) can be represented through quark distributions in the longitudinally — along 
the beam polarized nucleon 

9i{x)= E - > (3) 



For more details see 
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where qi+{x) and qi-{x) are distribution of quarks in the nucleon, polarized along or 
opposite to the nucleon spin, - are the quark charges. In the parton model x has 
the meaning of the part of the nucleon momentum, carried by quark in the infinite 
momentum frame (or for fast moving nucleon, e.g. in the virtual photon-nucleon 
c.m.s. at oo). Therefore, the integrals 



have simple physical meanings: they are equal to the parts of nucleon spin projection, 
in longitudinally polarized nucleon carried by the quarks of flavour i — u,d, s, ... . 

The structure functions Gi,G2 can be expressed in terms of the virtual photon 
absorption cross section (T3/2 and (T1/2 as well as the quantity o"/ describing the 
transition from transverse to longitudinal virtual photon polarization (on vice versa) 
in the forward scattering amplitude. The cross sections (73/2 and (T1/2 correspond 
to the projections 3/2 and 1/2 of the total photon nucleon spin upon the photon 
momentum direction. 

-V--3/. = ^^^(.G, + ,^G2) 

= 4 A Q/^ ^ j^^ X K ) 

In turn, (T3/2, (^1/2 and o"/ arc proportional to the s-channel helicity amplitudes 
Im Ti _i/2;i,-i/25 I'm Ti i/2;i,i/2 and Im Ti i/2;o,-i/2, where the flrst two indeces 
correspond to final photon and nucleon helicities, and the second two to initial ones. 
The interference cross section o"/ satisfies the inequality 

1(7/1 < ^JotOq = "/Rcto , (6) 

where ut and ctq are the absorption cross sections of transverse and longitudinal 
virtual photons in nonpolarized scattering, Gt = {0-3/2 + o'i/2)/2, R = (Jq/(Jt- 

The calculations of the structure functions in QCD are performed in the frame- 
work of the operator product expansion (OPE). The operators are classified accord- 
ing to twist = dimension - spin. The expansion in twists is equivalent to expansion 
in inverse powers of at fixed x. This expansion starts with twist 2 terms in the 
case of Gi, and with twist 3 terms in the case of G2{x, Q^), corresponding of scaling 
behavior of gi{x,Q'^), g2{x,Q'^). In perturbative QCD only the evolution with 
of the structure functions can be determined, the values of the structure functions 
at some fixed = are taken from the experiment and used as an input for 
evolution equations. 

The goal of all performed till now experiments was to measure the structure 
function gi{x, Q^). This goal was achieved in the following way. The experimentally 
measurable quantity is the asymmetry 

A - (7) 
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where a^^ , a^^ correspond to the cases, when the spins of longitudinal polarized /x(e) 
and nucleon are antiparallel or parallel. The asymmetry A is related to 



A = ((7i/2 - (73/2)/((Ji/2 + (J3/2) , A2 = (Ji/cTt (8) 



by 

where 



A = D{A, + riA2) (9) 

D= Q2(l-^ 

y2 + 2il-y)il + R) ' ' E y{2-y) ^ ^ 

y = E — E', E and are in initial and final fi{e) energies. It can be shown using 
the inequality (6), that the second term in (9) is small in practically interesting 
cases. Then the measurable quantity is proportional to Ai and the latter determines 
gi{x, Q^) through the relation 

where ^2(2;, Q"^) is the structure function for nonpolarized DIS. Before going to new 
developments in the theory of polarized structure functions, I recall some facts which 
should be well known. However, sometimes a misunderstanding occurs even here. 

3. Small X Domain — Regge Behavior 

Consider the spin-dependent virtual photon-nucleon forward scattering ampli- 
tude 

T'^uiP^ Q) = '^^nuXaqx !^^SaSi{i^, q^) + [i^s^ - {sq)pa] S2{iy, , (12) 

Gi(u,q'') = ImSi,2(i^,q^) . (13) 
The functions Si^2{i^,q^) satisfy the crossing-symmetry relations 

S,{v,q^) = S,{-iy,q'), S2{iy,q') = - S2{-iy,q') . (14) 

It can be shown, that at u ^ 00 and = Const (or x ^ 0,Q^ — Const) Si is 
contributed by Regge poles, satisfying the condition 

^(-l)^^ = - 1 , (15) 

where G is G-parity, T - is the isospin and a is the signature. Then in accord with 
(12), (14) a — —1 and the lowest Regge trajectory contributing to Si are 01(^41) 
and fi{D) meson trajectories 

Si{i^,q') ~ PaM^'""'-' , 9i{x) ~ . (16) 
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The intercept «„! is known not quite well, certainly it is negative and about ~ 
— (0.0 — 0.3). The same intercept is expected for /i trajectory. 

Up to the small contribution <iy~^'^, arising from (see eq.(5)), the function 5*2 
is proportional to the helicity amplitude Ti i/2;o,-i/2- The Regge poles do not con- 
tribute to the latter. This statement directly follows from the factorization theorem 
for Regge poles. Indeed, the interaction vertex of virtual 7 with any Regge pole in 
this case is a scalar, constructed from vectors q, e^, e^, hnear in the last two, what 
is impossible. For this reason 82(1', q^) at large u is determined by branch points 
contributions (three pomeron cut etc.) 

S2{i^,q')^Pc{Q')/ln'u+ E P^{Q>'''-'/lnu + ... . (17) 

i=P',A2 



4. The Parton Model 

As was already mentioned, the spin dependent structure function gi{x) is rep- 
resented in the parton model in terms of quark distributions by eq.3 in comparison 
with the representation of nonpolarized structure function 

F2{x)=xJ2e^[qi+{x)+qi4x)] . (18) 

i 

The standard paramctrization of quark distributions at some normalization point 
— Qo (usually Qo ~ 5 GeV"^) used in the descripton on the nonpolazed DIS data 
is obtained by matching the behavior at small x (Regge domain) with the behavior 
at large x, 1 — x <^ 1, following from the quark counting rules: 

q{x) = q+{x) + q.{x) = Ax~'^{l - xf . (19) 

It is expected that for valence quarks a ^ 0.5 (p - intercept), f3 ~ 3.0, for sea quarks 
a ^ 1.0 (experimentally a ~ 1.2 is more preferable), /3 ~ 5. 

A more refined parametrization, which accounts for the fact that q+{x) and q-{x) 
behave differently at x — > 1 q-{x)/q+{x) ~ {1 — xY was suggested by Brodsky (see 
the review ^). 

The same parametrization (19) can be applied to spin dependent quark distri- 
butions 

q+{x)-q_{x) = Bx-^'il-xf , (20) 

where a' = ~ —(0.0 — 0.3), ~ 3 for valence and ~ 5 for sea quarks. 

The parametrizations (19), (20) results in an interesting inequality in the case 
of strange (sea) quarks For strange quarks (19), (20) reduce to 

Xs{x) = x[s+{x) + S_{x)] ^ ^^(1 - x)^" , .^j^N 

s+{x) — s_{x) ^ Bs{l - x)/^" , ^ ' 
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where it was accepted for simplicity that a = 1.0, a' = 0. Since s±{x) > 0, we have 
an inequaUty 

As > \Bs\ , (22) 

which results in 

1 1 

V2,s = j xs{x)dx > J dx[s+{x) — S-{x)] = As (23) 



in the notation (4). The inequality (23) states, that, if the parametrization (23) is 
correct at some Ql, then at this Ql the part of proton spin projection carried by 
strange quarks cannot exceed the part of proton momentum carried by them. It 
is clear that small deviation from a = 1.0, a' = 0.0, say Acty = —0.2, Aa' — 0.3 
cannot seriously violate the inequality (23). 

5. Sum Rules in the Pcirton Model 

The sum rules for spin dependent structure functions of deep inelastic scattering 
follows from space-time representation of amplitudes in terms of current commuta- 
tors: 

1 

4 

Current commutators satisfy the causality condition 



T^^X) = ^ < P, s\\j^{x),j^{Q)\\p, S >antisym. ■ (24) 



b-^(x),>(0)] = at x^<0. (25) 

It can be shown that the domain near light cone corresponds to large Q'^ ^. The sum 
rules can be derived by considering cq.(24) in the limit x'^ ^ 0,Xq ^ and using 
equal time commutation relations for calculation of current commutators, expressed 
in terms of free quark fields. Accounting for m, d, s quarks we have from consideration 
of space-like components fi = i^v — k m. (24): 



Im Tjfc(x)^o^o = -^iki <p,s\- 



+ IjIi \p,s> , (26) 



where jii,j^i and j^i are isovector, octet and singlet (in SU(3) flavour symmetry) 
axial currents. The matrix elements in the r.h.s. of (26) are proportional to the 
nucleon coupling constants with axial currents. The isospin symmetry determines 
the proton matrix element of isovector current 

<P,s\ i|(0) \p,s>p = -2mgASi , (27) 

where Qa = 1.25 is the neutron /?-dccay axial coupling constant. On the other side, 
using the Fourier transformation of cq.l, the l.h.s. of (26) can be represented in 
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terms of the structure functions (for details see In this way the famous Bjorken 
sum rule ^ arises 

} 1 

j [gf{x) - g-{x)]dx = -gA . (28) 



The Ellis- Jaffe sum rule ^ for octet current follows if SU(3) flavour symmetry for 
baryonic octet /3-decays is supposed 

J\l{x)dx = ^(3F-L>) , (29) 

where F and D arc /3-decay axial coupling constants in the baryonic octet. 

In a similar way, by considering the behaviour of /mT^^(a;) (24) near the tip of 
light cone the Burkhardt-Cottingham ^° sum rule for the structure function g2{x) 
can be derived ^ 

/ g2(x)dx = . (30) 
Jo 

In the derivation (30) no current algebra is used, only the causality condition (25) 
and the hypothesis of absence of nontractable singularities in the amplitude at the 
tip of the light cone are imposed. 

The sum rule (30) looks like a super convergent sum rule, as it was originally 
derived in Ref.^*^. This is not the case, however. Indeed, in accord with (17) we can 
write the subtractionless dispersion relation for 5*2(2^, q"^) 

Q2/2 

If 82(1^, q^) would decrease faster than 1/z/ at i/ — > 00, then from (31) the sum rule 
would follow 

00 

j dv'G2{v\q^) = 0, (32) 

which goes into (30) in the limit Q'^ ^ 00. But due to contributions of pomeron, P' 
and A2 cuts, 82(1', q"^) decreases slower than l/u (see (17)) and (32) is invalid. That 
means, that there are higher twist contributions to g2{x, Q^) strongly increasing at 
small X in comparison with lowest twist one, of the order {l/Q'^){l/x'^)ln~'"x, n>5. 
(The sum rule (30) validity was confirmed by QCD sum rule calculation in the lowest 
order of the operator product expansion in - nucleon virtuality ^^.) 

Finally, I discuss the Gerasimov, Drell, Hearn (GDH) sum rule GDH sum 
rule refers to the scattering of the real photon on polarized nucleon. But, as we will 
see, it is closely related to the sum rules for DIS. 
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In the case of real photon only one amplitude 5*1 (z/, 0) survives in (12), because 
the kinematical factor in front of the other - 5*2 (z/, 0) is zero. The function S'i(z/, 0) 
satisfies the unsubtracted dispersion relation 



oo 



S,(u,0) = Afv'dv'^lp^^. (33) 



Let us go to the limit z/ ^ in this relation. According to the F.Low theorem, 
the terms, proportional to zero and first power in the photon frequency in the 
photon-nucleon scattering amplitude are expressed via the static characteristics of 
the nucleon. The direct calculation of Feynman graphs gives 

^i(0,0)p,„ = (34) 

where ^ are the anomalous proton and neutron magnetic moments. Prom (33), 
(34) the GDH sum rule follows 



°o , 1 
-^Gi,,,„(z.,0) = --<„. (35) 



/ 





The essential point is that at = and z/ ^ the spin dependent forward 
Compton amplitude is a constant and has no nucleon pole. That means that in 
the l.h.s. of the sum rule (35) the elastic contribution is absent — only inelastic 
processes contribute. That is why the GDH sum rule is very nontrivial. 

6. QCD Corrections to Sum Rules 

In what follows we use the notation 

1 



In the parton model 



1 / 1 4 

Tp.n = — (^iaa + -as + -S ) , (37) 



where, for 3 flavours in notation of (4) 



= Au- Ad ; as = Au + Ad - 2 As 
E = Alt + Ad + As 



(38) 



and E has the meaning of part of the nucleon spin projection carried by quarks. 
The sum rules, presented in Sec.4 in the supposition of SU(3) flavour symmetry in 
baryonic /3-decays result in equations 

as ^ 9a , as = 3F - D . (39) 
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In QCD eq.37 is modified by account of as-corrections. Tlie first order corrections 
in as were calculated in i3-i7^ ^j^g ^,2 ^,3 corrections to the octet part of the 
sum rule in the corrections to the singlet part in With QCD corrections, 
calculated in MS regularization scheme the sum rules takes the form (for 3 flavours) 



1 

12 



1 ^^-^ - 3.6 



TT 



TT 



X 



(ias + \a,) + I 



3 



TT 



20 



TT 



TT 



X 



S - 



where Nf — 3 is the number of flavours, 

Ag{Q')^ f dx[g^{x) - g^{x)] , 
Jo 



(40) 



(41) 



5f±(a;) are the distrubutions of gluons in the longitudinally polarized nucleon with 
spin projections along or opposite to the nucleon spin. /S.g has the meaning of the 
part of nucleon spin projection carried by gluons. In eq.(40) 03,03 and E are 
independent. 

A physically new point in eq.40 is the appearance of gluonic contribution to the 
sum rule. Although the last term in (40) is an ag correction and naively it would 
be expected to vanish in the hmit Q"^ — > 00, in fact it does not, because the Ag{Q'^) 
anomalous dimension is equal to —1, i.e. 



and 



lim as{Q^)Ag{Q^) = Const 

>oo 



(42) 



(43) 



For this reason the last term in (40) cannot be considered as a correction vanishing 
at — > 00, like any other as{Q^) corrections (e.g. the ~ ccg terms in the figure 
bracket in (40), where 03,08 and E have zero anomalous dimensions). 

On the other side from the nucleon spin projection conservation it follows that 



1E + A^(Q2) + L,(g2) = l , 



(44) 



where Lz is the orbital momenta contribution. Therefore, the logarithmically in- 
creasing with Ag{Q'^) means that L^iQ'^) ~ InQ"^ at large Q^, is negative and 
compensates Ag{Q'^). As a direct consequence the quark model of nucleon with 
quarks in 5— states cannot work at high Q^. 
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There was a wide discussion in the past years if gluonic contribution AgTp^n to 
Tp^n (the last term in eq.(40) ) is uniquelly defined theoretically or is not ^''''20-32 
The problem is that gluonic contribution to the structure functions, described by 
imaginary part of the forward '-fmrt - gluon scattering amplitude (Fig.l) is infrared 
dependent. Since in the infrared domain the gluonic and sea quark distributions 



Fig. 1 

The photon-gluon scattering 
diagrams, the wavy, dashed 
and solid lines correspond to 
virtual photons, gluons and 
quarks. 



are mixed and their separation depends on the infrared regularization scheme, a 
suspicion arises that A^Fp „ can have any value. This suspicion is supported by 
the fact that in the lowest order in the terms, proportional to InQ"^ are absent 
in AgTp^n and this contribution looks like next to leading terms in nonpolarized 
structure functions, where such an uncertainty is well known. 

In order to discuss the problem consider the gluonic contribution gip{x, Q^)gi to 
the proton structure function gip{x, Q^), described by evolution equation 

9i,{x, Q\i = Nf^^ f (^] [g^{y) - g_{y)] , (45) 

where the asymmetry A{xi) is determined by the diagrams of Fig.l. The calculation 
of the asymmetry A{xi),Xi — —q^/2pq results in the appearance of integrals 

f (fk± 

1 ~7 \~n 5 ^ — -^7 ^ 7 

J (k'j_ — xi{l — xi)p'^ + rrig] 

which are infrared dependent. To overcome this problem it is necessary to introduce 
the infrared cut-off (or infrared regularization), to separate the domain of large k\, 
where perturbative QCD is reliable, from the domain of small k'j_. The contribution 
of the latter must be addressed to noncalculable in perturbative QCD parton distri- 
bution. Such a procedure is legitimate because of factorization theorem, which states 
that the virtual photoabsorption cross section on the hadronic target /i, (7^(x, Q^) 
can be written down in the convolution form 

Q') = E ^7(^, Q', M') ® /,/,(x, M') , (46) 
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where a] is the photoproduction cross section on i*^ parton {i = q,q,g), fi/h sue 
the parton distributions in a hadron h, (S> stands for convolution. Both a] and 
fi/h depend on the infrared cut-off M^, but the physical cross section is cut-off 
independent. The variation of corresponds to redistribution among partons: the 
trade of gluon for sea quarks. 
As follows from (45) 



< > 



A.r^ = Nf^^AgA{M^) , (47) 

where ^ 

A{M^= f dxiA(xuM^) . (48) 
Jo 

The convenient way is to introduce cut-off in 

kl > M2(xi,p2) (49) 

Generally, may depend on xi and p^. For example, the cut-off in quark virtuality 
in Fig.l -/c^ > corresponds to the form of (49) with = (1 - Xi){M^ +P^Xi) 
if xi < -Ml/p^. 

In the calculation of the diagrams of Fig.l it is reasonable to neglect the light 
quark masses in comparison with gluon virtuality p^, since we expect that | | is 
of order of characteristic hadronic masses, |~ 1 GeV'^ Then introducing the 
infrared cut-off (49) we have 



A(M2) = |l -y'cia;i(l -2xi)[Znr-r]| , 



(50) 



where 

xi(l-xi)p-M2(xi,p2) ■ ^ > 

If = Const - a rectangular cut-off in A;^^ - , the integral in (50) vanishes (the 
integrand is antisymmetric under substitution xi — > 1—xi). Then A — —as/^i: and 
we obtain the gluonic contribution to Fp^^ (40). However, other forms of M'^{xi,p'^) 
result in different values of A(M^), what support the claim "^-^ that A(M^) is cut-off 
dependent. Even more, if we put = 0,m^ 7^ - the standard regularization 
scheme in the calculation of nonpolarized deep inelastic scattering - we will find 
A — O"^^. This result is, however, nonphysical because the compensation of —as/27r 
term in A arises from soft non-perturbative domain of k\ ~ m^, which must be 
attributed to sea quark distribution. 

Although generally A, as well as Ag and AgTp, are infrared cut-off dependent, 
a special class of preferable cut-off's can be chosen. In OPE the mean asymmetry 
A is proportional to the one-gluon matrix element of axial current 

^nXa{0:P:P) = < fi-l^A | jn5{0) \ gi£„ > , (52) 



10 



at zero momentum transfer. It can be shown that this quantity is proportional 
to the divergency lijT^i\a{i,Pi,P2)J = P2 — Pi, in the hmit P — > 0, i.e. to the one 
gluon matrix element of axial anomaly. The standard expression for axial anomaly 
corresponds to = Const and A — —as/^Ti (although it was demonstrated 
that the same infrared dependence persists here also). Therefore, the cut-off in < 
= Const is preferable, since it preserves the standard form of axial anomaly. In 
this case the axial anomaly can be considered as a local probe of gluon helicity. Some 
care, however, is necessary, when the calculations with this cut-off are compared with 
HO calculation in nonpolarized scattering, where, as a rule, another regularization 
procedure — = 0, 7^ - is used. 

The gluonic contribution to the spin dependent proton structure function 
gip{x, Q^)gi can be found experimentally by measuring inclusive two jets produc- 
tion in polarized DIS However, only large - k\ component of Fig.l diagrams 
can be determined in this way, the small - component is unmcasurable, since in 
this case it is impossible to separate one jet from two jets events This fact is 
in complete accord with the formulated above statement about the arbitrariness of 
infrared cut-off. 

What numerical value of gluonic contribution to Tp^n{Q'^) can be expected? 
Naively, we can say, that since gluons are carrying about 50% of proton momen- 
tum < X >gi— 0.5 * , they can carry the same amount of proton spin projec- 
tion, Ag ^ 1/4 — 1/2 (see (44) ). The same conclusion follows from the simplest 
parametrization of gluon distribution like (19), (20), or from the parametrization, 
suggested by Brodsky 1 If Ag ^ 1/2 at = 10 GeV^ (^^(lO GeV^) ^ 0.25 at 
Aqcd = 200Mey) then 

^g^p,n ~ 0.0066 . (53) 
Earlier Brodsky and Schmidt basing on the model of bound-state nucleon and 
counting rules suggested a different parametrization and obtained Ag ^ 1.2 at 
~ 1 GeV^. Up to terms of order the perturbative evolution of Ag{Q^) with 
Q"^ is given by ^^'^^ 

where b = 11 — {2/3)Nf = - is the normalization point. As follows from (54), 
if the value = 0.25 is accepted at /x^ = IGeV'^, then Ag{10GeV^) 0.6 (at 
^QCD = 200MeV and S(//^) = 0.59 - see Sec. 6). In order to have higher values of 
Ag (lOGeV^^) it is necessary to start from lower normalization point. 



This number refers to typical Q"^ in DIS experiments, ~ 5 GeV'^. At 
>2 GeV^ < X >gi only slighly depends on Q'^, but is decreasing at low 

up to < X >gi^ 0.3 at = 0.3 GeV'^ Since Ag{Q'^) increase as In at 
high Q^, the conclusion Ag < 0.5 is valid if we assume, that gluonic distribution is 
described by equations like (19), (20) at ^ 2 - 10 GeV^. 
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E.g. (WGeV^) = 3 requires Ag = 1 at = O.SGeV^ ^5 rj.^^ ^^^.gg ^^^^^ 
Ag in the latter case is achieved due to increasing of g+{x) — g-{x) at small x and 
the problem arises if such increasing is compatible with experiment and does it pro- 
ceeds in the domain of x, measured in the existing experiments. There are different 
opinions about this subject 30,3i,35-37 pg^-j^aps the problem must be reinvestigated 
in the light of new experimental data. 

It is interesting to mention that in the ratio of Gross-Llewellyn-Smith and 
Bjorken sum rules the perturbative QCD corrections cancel up to (see also 

the discTission of this point in ^). Therefore the experimental examination of this 
prediction at moderate Q'^ could be also a good check of perturbative QCD and the 
role of higher twist terms. 

Let us turn now to twist four corrections to the sum rules (40). The general 
theory of these corrections was developed by Shuryak and Vainstein According 
to this theory the twist 4 corrections to (40) are proportional to the one-nucleon 
matrix elements of the operators 

< N I t/^^'^^ I AT >=< TV I gqGl^^^h'^il, Ts)q \N> , (55) 

where (5^^ is the dual gluonic field tensor, S{NS) are singlet (nonsinglet) in flavour 
operators, 1(t3) in the r.h.s. of (55) correspond to S{NS). In (55) q means u and 
d and the contribution of strange quarks is omitted. An attempt to determine the 
matrix elements (55) in the framework of QCD sum rule 3-point function calculations 
was done by Bahtsky, Braun, Kolesnichenko Their final results are 

1 ^,0.09 ±0.06^ 2m2 /• ^ , 

AtwisuTp,n = -3^(4, 1) + J ^ 9l,P,n{^)d^ ■ (56) 



The magnitude of twist 4 correction following from (56) is very small. Even with 
account of (negative) error it comprises ~ —1 • 10^^ at Q'^ — lOGeV'^ for proton as 
well as for neutron. 

Unfortunately, I cannot believe that the result (56) is reliable for the following 
reasons. In finding the vacuum expectation values induced by external axial field 
- a quantity, which essentially determines the final answer - the authors of ref. 
have really taken the octet field instead of the singlet one and used the dominance 
of massless goldstones (vr or 77) which is incorrect for the singlet field case. For 
the nonsinglet case the situation is, in principle better. But, unfortunately, in this 
case the main contribution, to the sum rule used by the authors, comes from the 
highest order accounted term in the OPE - the term ~ mg < t/jt/j >^ of dimension 
8. It is possible that in this case the contribution of unaccounted terms of OPE is 
essential or even that the OPE series is divergent at characteristic values of the Borel 
parameter = 1 GeV"^ used in the sum rule. On the physical side of the sum rule 
the contribution of continuum comprises 80% and the nucleon pole term from which 
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Tp^n was obtained gives only 20%. This also spoils the accuracy of the calculations. 
Finally, the result of (unlike the other results obtained by the QCD sum rule) 
depends on the ultraviolet cut-off. This circumstance introduces a noncontroUable 
uncertainty into the calculation. 

7. Nucleon Axial Coupling Constants 

As was demonstrated in the previous Sections, according to the sum rules the 
integrals 

are expressed through the nucleon matrix elements of isovector, octet and singlet 
axial currents, 03, ag and S. The magnitude of 03, determining Bjorken sum rule, is 
well known. According to isospin invariance it is equal to the axial coupling constant 
in neutron /3- decay, = = 1.257 ± 0.003 The theoretical accuracy in this 
number is given by the accuracy of isospin symmetry, i.e. it is of order 1%. In a 
similar way as is proportional to the matrix element of octet axial current 

ag = <P,s\ j^5^ \p,s> /2m . (57) 

This matrix element can be found if SU{3) flavour symmetry in baryonic octet 
/9-decays is assumed. Then 

ag = 3F - D = 0.59 ± 0.02 (58) 

and the numerical value in (58) corresponds to the best fit of neutron and hyperon 
/^-decays in SU{3) symmetry. The use of SU(3) symmetry in determination of og 
was questioned by Lipkin Indeed, the nondiagonal matrix elements enter neutron 
and hyperon (3 -decays, while the diagonal one appears in (57). If SU{3) symmetry 
is badly violated, it is possible that the latter has nothing to do with the formers. 
The accuracy of SU (3) symmetry in the matrix elements of octet axial current can 
be estimated in the following way ^. In SU{3) symmetry the interesting for us 
combination 3F — D can be found from any pair of /3-decays. If SU (3) is strongly 
violated we would expect that the values, determined from various pairs of /3-decays 
are different. Experiments tell us 

{3F - D)np,Ap = -3gA,np + Gqaap = 0.537 ± 0.09 , 

(3F - D)np,j:n = 9A,np + 2^A,Sn = 0.577 ± 0.034 , 
(3F - L')Ap,En = {^/2){gA,np + (7A,En) = 0.567 ± 0.034 , 

(3F - D)sA = 3gA,sA = 0.75 ± 0.15 . 

We see that all values (59) argee very well with (58) and, at least, the value of 
3F — D less than 0.5 is completely excluded. 
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The matrix element as (57) was also determined theoretically by QCD sum rule 
method In this calculation no SU (3) flavour symmetry of baryonic /^-decays was 
assumed. The result 

as = 0.5 ± 0.2 (60) 

is in agreement with (58). 

To estimate the proton matrix element of singlet axial current (or E) is a more 
complicated, till now nonsolved problem. Naively, by analogy with a part of the 
proton momentum, carried by quarks, which is about 50% at ~ 5 — 10 GeV^, 
we expect, that E 0.5. The same estimate follows also from the assumption that 
strange sea is nonpolarized (the original EUis-Jaffe ^ assumption). Then E f» 
given by (58). 

On the other side, assuming that the Skyrme model is suitable for the description 
of proton constituents, Brodsky, Ellis and Karliner argue that E is small of order 
1/Nc in the limit of large number of colours Nc- Indeed, in the Skyrme model 
the nucleon is built from octet of goldstones, which is decoupled from singlet axial 
current at Nc oo. The Skyrme model, probably, gives a good description of 
nucleon periphery. However, the description with this model of deep inside nucleon 
structure, which is measured in deep inelastic scattering, is questionable (see also 

In ref.48 an attempt was performed to calculate the nucleon matrix element of 
singlet axial current Oq = E by QCD sum rules in the similar way as it was done 
in the calculation of ag in The main difference with Og calculation is that in 
this case the anomaly plays an important role in the vacuum expectation value 
induced by external axial field, which essentially determines the final result. This 
problem was overcome in but nevertheless the result was negative: an attemt to 
determine E fails, since it was found that the OPE breaks down for singlet axial 
current in longitudinal and transverse parts of polarization operators and/or in the 
NNA vertex function. The physical consequence of this consideration is that one 
should expect o noticeable Okubo-Zweig-Iizuka (OZl) rule violation in the nonet of 
axial l"^"*" mesons. 

The authors of Refs. 24,49 tried to use the so called "C/(l) Goldberger- 
Treiman" relation 

Fno{0)grn,NN = '2mao (61) 

to determine oq- Unfortunately, no definite conclusion could be obtained in this 
approach without serious additinal hypothesis. The reason for this comes from the 
fact that due to anomaly rjo in (61) is not a goldstone and not a physical rj' meson. 
Perhaps, the troubles appearing in this approach are of the same origin as in 
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8. The Experimental Data on gi{x) and Their Interpretation 

The "Sturm und Drang" period in the investigation of nucleon spin structure 
started after appearance of famous EMC results ^, where the structure function 
gip{x) was measured in the interval 0.015 < x < 0.47 at the mean = 10.7 GeV^. 
Combining their data with earlier SLAC ^ data, where the measurements were made 
at 0.1 < X < 0.65 at ~ 3 GeV"^ and performing the extrapolation in the domain 
X < 0.015 and x > 0.65, EMC obtained 

Fp = 0.126 ± OmO{statist) ± 0.015 (system.) , = 10.7 GeV^ (62) 

Compare (62) with the simplest theoretical version — the Ellis-Jaffe sum rule, 
— eq.40, where the gluonic contribution is neglected and it is supposed that strange 
sea is nonpolarized. As = 0. Then E = and using (58), we have at as{Q^ — 
10.7 GeV^) = 0.25 {Aqcd = 200 MeV) 

T^-^ = 0.171 ± 0.004 , (63) 

where the error comes from the error in (58) and from the uncertainty in the 
correction, which could be of the order of term (besides of uncertainty in Aqcd)- 
The Ellis-Jaffe prediction disagrees with (62) by 2.6 standard deviations. 

Reject now the hypothesis As = 0, but still neglect gluonic contribution. Using 
experimental value (62) we can determine from (39), (40), (58) the values Aw, Ad, As 
and E. 

Aii = 0.785 ± 0.06 , Ad = -0.47 ±0.06 , As = -0.14 ± 0.06 , , . 

E = 0.17 ±0.17. 

The results (64) as originally was stressed in ^'^^ are very strange and unexpected: 
the part of proton spin projection carried by quarks is small and compatible with 
zero, strange quarks are carrying a large amount of proton spin in contradiction with 
the quark model and inequality (23). Indeed the experimental data on DIS give ^° 

V2,s = 0.026 ± 0.006 (65) 

and (23) is strongly violated. 

The account of gluons does not improve the situation essentially, if Ag pa 0.5 at 
Q"^ — 10 GeV'^. Using (53), we find that in this case As and E increase in comparison 
with (64) correspondingly by 0.02 and 0.07. If, we assume that Ag'(10 GeV'^) ^ 2, 
then 

As = -0.06 ±0.06 ; E = 0.42 ±0.17, (66) 

the inequality (23) is fullfilled within the errors and S does not contradict the naive 
expectations. The value A5f(10 GeV"^) = 2 can be achieved, if Ag{l GeV"^) = 1.2 
what was suggested in This solution means, however, that the constituent quark 
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model would be strongly violated at ~ 1 GeV^, since, according to (44) we would 
have L,{1 GeV^) ^ -1. 

Leaving aside the discussion of the scenario with large Ag, let us consider if it is 
possible to have large | As |, violating the inequality (23). Qualitatively, the problem 
can be formulated as follows: could it be, that in one experiment we see a small 
number of strange quark pairs in the proton and a much larger one in the other? The 
answer is since the number of ss pairs is not conserved, the question how many 
strange quark pairs are in the proton is not correctly formulated until one specify the 
operator, with which the proton strange content is measured. The measurements 
of a part of proton momentum, carried by s-quarks corresponds to the matrix 
element of the s-quark energy momentum tensor operator, ~< p \ 0^^^, \ p >, 
the measurement of proton spin projection, carried by s-quark corresponds to the 
s-quark axial current operator As ~< p | j^s \ P >■ The small value of and 
a large of As would mean, that in the case of the operator the transitions 
ss — ^ uu + (id are suppressed, while this is not the case for the j^^ operator — the 
OZI rule works well in the first case and is violated in the second. This circumstance 
is not surprising. In the meson nonets the situation is similar: the OZI rule works 
well in the case of vector and tensor mesons and is 100% violated in pseudoscalar 
nonet. The latter phenomenon was explained theoretically by attributing ss — > 
uu + dd transitions in mesons to nonperturbative QCD effects — the instantons 
In the field of instanton the transitions ss ^ uu + dd proceed in pseudoscalar and 
longitudinal axial channels, but are forbidden in vector or tensor channels. 

It is plausible that the same instanton mechanism results in large ss — > uu+dd 
transition in < p | j^^ \ P >, but not in < p | ^^^^ \ P >■ If it is indeed the case, then 
the two-component form of s-quark distribution is expected 

s+ix) + s_(x) = {Ai/x){l - x)^' + {A2/x){l - xY' , . . 

s+{x) - s^x) ^ B{1 - xf , ^ ' 

where /3s ^ and p' p 10 — 12. It is easy to see that for distributions (67) the 
inequality (23) is weakened to | As | < 8V2S and (64), (65) are not in contradiction. 

In the last year there was an impressive progress in experiment in the problem in 
view. SMC presented the results of new measurements of gip{x). The region of x, 
where the measurement were performed was extended to 0.003 < a; < 0.7. The SMC 
data are in a good agreement with previous EMC and SLAC data in the domain of 
X common for all experiments. However, at x < 0.015 the SMC points lay higher, 
than the EMC extrapolation. The SMC data are presented in Fig.2. SMC found 
for Tp at g2 = 10 GeV^ 

Fp = 0.136 ± 0.011 ± 0.011 . (68) 
In the combined analysis of all data - SLAC, EMC and SMC, it was obtained 

Tp = 0.142 ± 0.008 ± 0.011 . (68') 



16 



Fig. 2 

The SMC data in comparison with the the- 
ory. The sohd circles (right-hand axes) show 
xgip{x), the open boxes (left-hand axes) show 
Ix gip{x)dx. Only statistical errors are shown. 
The solid square shows the SMC result for Tp 
with statistical and systematic error combined 
in quadrature. The solid and dashed lines are 
theoretical predictions ^ and correspondingly 
(see text). The Ellis-Jaffe prediction (63) is 
marked by the cross. 
If gluonic contribution to Tp is neglected, then it follows: 
from (68): 

E = 0.26 ± 0.14 , As = - 0.12 ± 0.05 , (69) 

from (68') 

E = 0.31 ± 0.13 , As = - 0.10 ± 0.05 . (70) 

Quite recently the E143 group at SLAC presented the preliminary result of gip{x) 
measurements at 0.029 < x < 0.8 and = 3 GeV'^. E143 found 

Tp = 0.133 ± 0.004 ±0.012, ^ 3 GeV^ . (71) 

(In obtaining (71) E143 group introduced a correction for x < 0.029 region contri- 
bution by extrapolation of their data. This extrapolation curve however, lays below 
the direct SMC data points.) E143 data (71) result in 

E = 0.32 ± 0.12 , As = - 0.10 ± 0.04 . (72) 

The neutron structure function gin{x) was measured in two experiments: by SMC 
at CERN in polarized muon scattering on polarized deuterium and by E142 group 
at SLAC in polarized electron scattering on polarized ^He The mean value of 
0^ in the SMC experiment was Q'^ = 4.6 GeV^, in E142 experiment = 2 GeV^. 
After correction for D-wsive admixture in the deuteron, SMC found 

+ r„ = 0.050 ± 0.044 ± 0.033 , = 4.6 GeV"^ . (73) 

In order to determine r„ from (73) and Tp measurements the data must be taken 
at one common value — Qq- We will do such recalculation, in the same way as 
in ^'^ - accounting for pcrturbativc corrections given in (40) and twist 4 correction 
(56) and neglecting Ag contribution. As Ql we choose Ql — 10.7 GeV'^ . We get 

(rp + r„) |q2 

(Tp + r„)4.6 GeV^ 



(74) 
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and the corresponding to (73) value at Qq is 

+ r„ = 0.053 ± 0.046 ± 0.035 , Ql = 10.7 GeV'^ . (75) 
If the result of combined analysis (68') is taken for Tp then 

r„ = -0.089 ±0.04 ±0.04 (76) 

and 

Tp-Tn = 0.231 ± 0.045 ± 0.045 , Q"^ = Ql . (77) 

This experimental number may be compared with the Bjorken sum rule, calculated 
with account of perturbative QCD corrections up to 

(rp-r„)^^^^^ = 0.186 ±0.003, = g^ (78) 

The two numbers agree within the experimental errors, which, unfortunately, are 
rather large. 

The E142 experiment has the advantage that the spins of two protons in ^He 
are compensated and the experiment (up to small correction) gives directly gin{x). 
The experiment was done at 0.03 < a; < 0.6. Performing the extrapolation in the 
regions of small and large x the E142 group obtained 

= - 0.022 ± 0.011 , = ^GeV"^ (79) 

Ellis and Karhner using E142 data, but different extrapolation obtained instead 
Tin = —0.028. The value Fi^ close to the latter was obtained by SMC when they 
use their data for extrapolation. Transferring the E142 group data to the common 
value = Ql, we have from (79) 

r„(g2) ^ - 0.024 ±0.012 (-0.030) . (80) 

(In the parenthesis the value, following from analysis ^^'^^ is given.) Prom (68') and 
(80) it follows for the Bjorken sum rule 

Pp-P„ = 0.166 ±0.021 (0.172), Q^ = Ql , (81) 

in agreement with theoretical value (78). The values of proton spin projection carried 
by quarks can be also determined from E142 data (79). The results are the following 

Au = 0.88 ± 0.05 , Ad = -0.37 ± 0.05 , E = 0.47 ± 0.12 

As = -0.04 ± 0.05 . ^ ' 

(The difference between (79) and the value given by ^^'^^ is included in the error.) 

The results of all experiments, represented in terms of the part of proton spin pro- 
jection, carried by ?i, d and s-quarks arc compatible, although some spread in the 
results is also seen. (Particularely, the difference between the EMC and E142 data.) 
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Different experiments were done at different and, even more, eacli bin in x 
corresponds to eacfi own Q^. Tfie problem arises, if tfie account of dependence 
could change the results. Experimentally, no dependence in the asymmetry A 
(7) was observed, but the accuracy is not good enough — not better than 10% at 
high and 20% at low . In ref.59 by solving the evolution equations the problem 
was studied, how the fact, that each bin in x is measured at each own Q^, can affect 
the final result. It was found that this effect is small, less than the errors in the 
present experiments. 

Till now when examining the experimental data the higher twist corrections were 
disregarded, or accounted basing on the calculations where they are very small, 
much less than experimental errors. However, as was explained in Sec. 5, the results 
of are not convincing. Now I will discuss another way of determination of ^p^tikQ^) 
Q"^ dependence, based on the connection of V^^JyQ'^) with GDH sum rule. 

Following introduce the functions 



The dependence of Ip^niQ"^)-, and the difference Ip{Q^) In{Q^) is plotted 




(85) 



(84) 



(83) 



Fig. 3 



The connection of GDH sum rule with sum 
rules at high Q'^ - the qualitative dependence 
of /p(g2),4(Q2) and/p(g2)-4(Q2) on Q\ 
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qualitatively in Fig. 3. The case of Ip is especially interesting. At large Ip{Q^) is 
positive, but Bit Q'^ — it is negative and large, /p(0) = —0.72. This means that in 
the region Q^<1 GeV'^, Tp{Q'^) varies strongly and changes the sign. Such a behavior 
cannot be achieved by smooth extrapolation of perturbative QCD corrections and 
is an indication of large nonperturbative effects. 

In ref.6 the vector dominance based model was suggested, which interpolates 
Ip,n{Q'^) at intermediate Q^, has the correct asymptotics (84) at large and satisfies 
the GDH sum rule (85) at = 0. In constructing such model it must be taken into 
account, that the contribution of baryonic resonances are important in ip,n(Q^) at 
low . Therefore, the model for Ip^niQ"^) has the form 



(86) 



where Ip'^niQ'^) the contribution of baryonic resonances, known from the analysis 
of pion electroproduction experiments (up to mass W — 1.8 GeV) I'p^niQ^) 
defined by 



-'P,n 



1 + -^ 



2m2r-„ 



(0) 



(87) 



p,n ' 



Fp^^ means the value Fp^„ 
to determine 



where /s^ is the vector meson mass, /x^ = 0.6 GeV^ 
at large with higher twist terms excluded. For our purposes 
the power corrections we may consider F^^^ as Fp^„ at ~ 10 GeV^. The 
second term in (87) in the framework of VDM corresponds to the case when in 
forward virtual 7 Compton scattering amplitude both 7 interact with the nucleon 
through vector mesons. The first term in (87) corresponds to the case, when one 
7 interacts through vector meson and the other one directly. It is easy to see that 
in the model the asymptotic behavior of /p,n(Q^) and GDH sum rule is fullfilled. 
From the experimental data on pion photoproduction, the known values of Hp^n and 
experimental data of F^^^^ it was found 



0.43 ±0.10 , 



0.0 



+0.3 
-1.2 



(89) 



The large error in c„ results from large uncertainty in F„ in E142 experiment. 

The resonance contribution is absent (or small) in the region of and x, where 
the existing experiments were performed. Then the power corrections came from 
Jp^ in (86) and can be calculated using (87) - (89). The results are 



Q\GeV^) 

'pas /'pexper 

pip 
pas I'^exper 



2 

1.44 
1.30 



3 

1.29 
1.20 



4.6 10.5 
1.19 1.08 
1.13 1.06 



(90) 
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The results for Tp^n, E and As arising after introduction of power corrections (90) 
are given in the Table I. (The SMC r„ data are not included, because of large 
uncertainty in E and As). 



Table I 





EMC(p) 


SMC(p) 


E143(p) 


E142(n) 


r 


0.137 ±0.018 


0.147 ±0.016 


0.172 ±0.013 


-0.029 ±0.015 


E 


0.28 ±0.17 


0.37 ±0.14 


0.68 ±0.12 


0.46 ±0.13 


As 


-0.10 ±0.06 


-0.08 ±0.05 


+0.02 ±0.04 


-0.04 ±0.05 



All the data after introduction of power corrections are consistent with E 
0.4 — 0.5 and As —0.02 — 0.04 — the values close to naive expectations — and 
no large gluonic term is needed. The contribution of gluons was disregarded in the 
calculation of Au, Ad, As and AE, the results of which arc presented in (69), (70), 
(72), (82) as well as in Table 1. As was mentioned above, if Ag{10 GeV'^) = 0.5 the 
account of gluons would increase E and As by 0.07 and 0.02 corresponding and 
twice as large if Ag{lQ GeV'^) = 1, what, may be, is also acceptable. 

Till now I discussed mostly the integrals of g\{x) over x, which are known better 
theoretically than the x-dependence of gi{x). There are few theoretical models for 
gi{x) but, unfortunately, the successful description of the data was achieved by the 
price of use some experimental information as input. Fig. 2 displays the theoretical 
curve for gip{x), given by Brodsky on the base of matching the Regge behavior at 
small X with quark counting at large x. An essential point was that the matching 
was done separately for q+{x) and q_(x){q — u,d,s). The values Am, Ad, As from 
EMC experiment arc used as input. The same figure shows also gip{x) found in the 
QCD sum rule approach In this calculation no experimental input is used: the 
only parameter, which enters the calculation is the fundamental QCD parameter - 
the quark condensate ag < \ 'ipip | Unfortunately, the domain of x, where 
the results are reliable (this domain was determined in the calculation) is rather 
narrow, 0.5 < x > 0.8. 

9. Suggestions for Future Experiments 

1. The experimental study of the sum rule dependence. As was discussed in 
the previous Section, the part of the deviation from the theoretical expectations, as 
well as the difference between the data, obtained at different Q^, may be attributed 
to high twist effects. Also, the role of dependent corrections is not negligible 
in the problem in view. For this reason the precise measurements of the depen- 
dence of the sum rules would be very desirable. Such measurements are expected in 
the near future at SLAC (see ^^'^^), CERN and DESY (HERMES CoUaboration). 
Especially interesting would be the direct experimental test of the GDH sum rule 
and the study of dependence of Ip,niQ^) (83) in the domain < Q^<2GeV^. Up 
to now the GDH sum rule was checked in indirect way, using the data on the pho- 
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toproduction of baryonic resonances. It was found that if the GDH sum rule 
is decomposed into isoscalar and isovector components, then the isovector-isovector 
sum rule is well satisfied, isoscalar-isoscalar contribution is small. However, in the 
isovector-isoscalar sum rule the contribution of high mass resonances is essential, 
what deteriorates the accuracy of the results. The study of Ip^niQ'^) dependence 
would be important for checking of various models. Particularely, in the model it 
is expected, that IpiQ"^) is crossing zero at Q'^ ^ 0.5 GeV'^. 

2. The search for Ag. The straightforward way to determine g+{x) — g-{x) is 
by mearuring two jets in the virtual photon fragmentation region in the polarized 
fi{e)- proton deep inelastic scattering Such events arise from 7 — — ^ qq collisions 
and may be separated from other two jet events by their known dependence 
When the jets are well separated in k"^, there are no infrared problems and the cross 
section is predicted uniquUy in terms of g+{x) — g-{x). 

For discussion of other possibilities see the review and references herein. 

3. The measurements of As. Besides of As determination in future precise po- 
larized deep inelastic ;u(e)-nucleon scattering. As can also be measured in elastic up 
scattering by separating axial formf actor contribution at zero momentum transfer 
The axial formfactor arises from Z-boson exchange and at = is proportional 
to the matrix element 

< P I u^y^^^u - d^i^-i^d - s'Ju'j^s I p >= -2s^m{Au - Ad - As) . (91) 

Such an experiment gives an additional equation for determination of Au, Ad, As. 
The existing data indicate on large As 

As = -0.15 ± 0.08 . (92) 

The determination of axial formfactor Ga{0) from the data must be done with some 
care. The experiment is performed ai ^ and the extrapolation to = 
proceeds, assuming the dipole formula for the formfactor. 

However, due to chirality conservation, (93) cannot be correct in the whole domain 
of Q^. This follows from the requirement that at large the axial formfactor 
must coincide with vector one Therefore a deviation from the dipole fit (93) is 
expected at intermediate Q"^. This expectation is in accord with the discussed in Sec. 8 
nonperturbative effects (instantons) in flavour singlet axial formfactor. Surely, the 
observation of the deviation from the dipole fit in singlet axial would be extremingly 
important for our understanding of nonperturbative QCD. 

The measurements of s-quark distribution in nonpolarized nucleon would be also 
helpful in understanding the origin of the strange sea. The s-quark distribution 
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obtained from the data ™ on charm production in uN scattering shows a more steep 
1 — X dependence 

xs{x) ~ (1 - x)^ , /3 6.5 , 

then the quark counting rule expectation (5 ^ h. This may be an indication in favour 
of nonperturbative (instanton) component (67) in s-quark distribution, but more 

precise dat needed. 

If As is large and negative, then a surplus of jets with leading strange particles 
can be expected in the direction opposite to the proton spin. This also could be 
tested. 

4. The check of the OZI rule in axial mesons. A connected problem is the de- 
termination if OZI rule holds in the nonet of axial mesons. If As is large and 
nonperturbative effects are important in singlet axial channel, then it is expected 
a noticeable OZI rule violation in the nonet of axial 1"*"^ mesons. There are 
indications, that it is indeed, what happened 

10. The Structure Function g2{x) 

Since there are very good reviews ^^'^^ on this subject and almost nothing new 
happened after their appearance, I will indicate here the basic points only. The 
structure function g2{x) can be measured in deep inelastic ^(e)-nuclcon scattering, 
when ii[e) is polarized longitudinally and nuclcon transverse to the beam direction. 
In this case the contributions of gi and g2 to spin asymmetry (measured, e.g. by 
changing the direction of nucleon spin to the opposite) are of the same order, but 
suppressed by the factor ~ 1/Q in comparison with the gi contribution in longi- 
tudinally polarized scattering. In the cases of other directions of nucleon spin g2 
contribution is smaller, than gi by a factor ^ 1/Q. This fact demonstrates that g2 
may be considered as a twist-3 structure function. 

For free particle g2 = 0. This is an indication, that g2 cannot be described in 
terms of parton distributions, because in the parton model we consider deep inelastic 
scattering, as proceeding on free partons. The other argument in the same direction 
can also be presented. Let us suppose, following Feynman that the transversely 
polarized nucleon in the infinite momentum frame is described by parton densities 
tf{x) with polarizations ±1 upon nucleon polarization direction. Then, it can be 
shown ^, that 

Here /x^, Cj and p\j are the mass, charge and mean square of the transverse momen- 
tum of parton j. At large Q, since are the current masses, the r.h.s. of (94) is 
small. But the integral over x of the l.h.s. is not small, according to sum rules, so a 
contradiction arises. The quark-gluon interaction and transverse momentum quark 
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distribution in nucleon are of importance in g2{x) determination. Both of them are 
beyond the parton model. For this reason our theoretical knowledge of g2 {x) is very 
scarce now. 

In OPE on the light cone g2{x) can be decomposed into two sets of operators 
The first is twist-2 operators, the same, which appear in gi{x) decomposition, 

the second is twist-3 operators. The first set contribution can be expressed through 

gi{x) and the relation can be obtained 

92{x, Q') = gY'^ix, Q') + Q') , g^'^ix, Q') = -g,{x, Q') + / -g^{z, Q') , 

J z 

X 

(95) 

where g^^ and g2 correspond to twist-2 and 3. There is no reason to believe, that 
the most interesting piece — g2{x,Q'^), determined by quark-gluon interaction, is 
much less than g'^^{x, Q"^). Therefore, it is unclear, if decomposition (95) is helpful. 

It is very important to check experimentally the Burkhardt-Cottingham (BC) 
sum rule (30). At such investigation it would be interesting to study dependence 
of the sum rule: as was argued in Sec. 4 one may expect that BC sum rule is valid 
only in low twist order and at ~ 1 — 2 GeV^ the violation of BC sum rule could 
be seen. 

11. Chirality Violating Structure Function hi{x) 

As is well known, all structure functions of the twist two - Fi{x), F2{x), gi{x) 
which are measured in the deep-inclastic Icpton-nuclcon scattering conserve chiral- 
ity. Ralston and Sopcr ''^ first demonstrated that besides these structure functions, 
there exists the twist-two chirality violating nucleon structure function hi{x). This 
structure function does not manifest itself in the deep inelastic lepton-hadron scat- 
tering, but can be measured in the Drell-Yan process with both beam and target 
transversally polarized. The reason of this circumstance is the following. The cross 
section of the deep inelastic electron(muon)-hadron scattering is proportional to the 
imaginary part of the forward virtual photon-hadron scattering amplitude. At high 
photon virtuality the quark Compton amplitude dominates, where the photon is 
absorbed and emitted by the same quark (Fig.4a) and the conservation of chiral- 
ity is evident. The cross section of the Drell-Yan process can be represented as an 
imaginary part of the diagram. Fig. 4b. Here virtual photons interact with different 
quarks and it is possible, as is shown in Fig. 4b, that chirality violating amplitude in 
Drell-Yan processes is not suppressed at high Q'^ in comparison with chirahty con- 
serving ones, and consequently, corresponds to twist two. However, this amplitude, 
corresponding to target spin flip, has no parton interpretation in terms of quark 
distritutions in the helicity basis and, as was shown by Jaffe and Ji can be only 
represented as an element of the quark-quark density matrix in this basis. 
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Fig. 4 

a) Deep inelastic lepton-hadron scatter- 
ing, the quark chiralities are conserved. 
Solid lines are quarks, wavy lines are vir- 
tual photons, R{L) denote right (left) chi- 
rality of quarks; b) Drell-Yan process with 
chirality of quarks flipped. 

However, hi{x) can be interpeted as a difference of quark densities with the 
eigenvalues -1-1/2 and —1/2 of the transverse Paili-Lubanski spin operator .§±75 in 
the transversely polarized proton. It this basis hi{x) can be described in terms of 
standard parton language. 

Until now there are no experimental data on the chirality violating nucleon struc- 
ture function hi{x). Besides ^^'^^ the theoretical study of this structure function have 
been also performed by Artru and Mckhfi The role of hi{x) in factorization of 
a general hard process with polarized particles was investigated by Collins The 
first attempt to calculate hi (x) was carried out by Jaffe and Ji by means of the 
bag model. 

The proton structure function hi{x) can be defined in the light cone formalism 
as follows 

^ jge^^^ <p,s I ^(0)(T^^75^/'(An) \ p, s >^ 2[hi{x,Q^){s±^p^ - s±^p^,)+ 
+hL{x, Q'^)m'^{pfin^ - p„nn){sri) + h^{x, Q'^)m?{s^ijn„ - s^_i,n^)\ . 

Here n is a light cone vector of dimension (mass)~^,n^ = Q,n^ = Q,pn = l,p 
and s are the proton momentum and spin vectors, = rn?,s^ — — l,ps = 
and s — {sn)p -\- {sp)n -\- s±,hL{x,Q^) and h3{x,Q^) are twist-3 and 4 structure 
functions. For comparison in the same light cone notation the standard structure 
function Fi{x,Q'^) is given by 

r dX 

J 2^"^'^" ^ ^' ^ ' V'(0)7^V^(An) \p,s>= 4[Fi{x, Q^)p^ + M^U{x, Q^)n^] (97) 

(Eq. (96), (97) are written for one flavour). 

Basing on the definition (96) an inequality was proved in 

q{x) > hl{x) , (98) 

which holds for each flavour q = u. d. s. (Here h^x) is the flavour q contribution to 
hi{x).) It was also suggested that | hi{x) \>\ gi{x) \. 

hi{x) can be also represented through a T-product of currents 

T,{p,q,s)^i J d'xe"^^ <p,s\{l/2)T{j,,{x),m+j{x),j,,{0)}\p,s> , (99) 
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where j^^^x) and j{x) are axial and scalar currents. 
The general form of ^^(p, s) is 



T,{p. q, s) = i^s, ^g^j h,{x, Q') + (^p, - ^ j (qsMx, Q') + 

+ £tJ,v\aPuq\Sa{qs)l2{x,q^) (100) 
(only spin-dependent terms are retained). It can be proved ™, that 

/ii(x,Q2) ^ --Imhx{x,Q'^). (101) 

TT 

As is clear from (96) or (99) hi{x) indeed violates chirality. It is worth mentioning 
that instead of axial and scalar currents in (99) the combination of vector and 
pseudoscalar currents jn{x)js{Q) — js{x)j^{0) can also be used (notice the minus sign 
in the crossing term). 

The structure function hi{x) can be measured in the production of jJi^ ^~ or 
e+e~ pairs in pp (or pp) collisions, when both protons (or p and p) are polarized 
transversally to the beam direction. Introduce the coordinate system, where the 
beam direction is along z, the proton spin direction — along x and the momentum 
of the lepton pair is characterized by polar and azimuthal angles 6 and </? (in the 
lepton pair c.m.s.). Then the asymmetry is given by 

^ a+{x,y) - a^{x,y) ^ sin^cos2(^ Eg eX(^)^f (l/) r-^^^] 
^ a+{x,y)-a^{x,y) {I + cos' 9) a elq^{x)qi{y) ^ > 

Here (T+(a;, y) and (T_(a;, y) are the cross sections of lepton pair production, where 
the spin of one of protons is along or opposite to spin of the other (or p) . 

In (102) X and y are the part of the momenta, carried by quarks and antiquarks, 
annihilating into virtual photon with momentum q 



q = xpi + yp2 + q± 
M'^ — q' = xys — q\ 



2 _ „2 _ .2 (103) 



where M is the mass of leptonic pair, and are the structure functions hi for 
quarks and antiquarks in the target and projectile correspondingly, q"- and q°' are 
the quark and antiquark distributions. In pp collisions h\ = hl^q"' — q"' . 

Since hi{x) violates chirality one may expect that in QCD it can be expressed 
in terms of chirality violating fundamental parameters of the theory, the simplest 
of which (of lowest dimension) is the quark condensate. Basing on this idea hi{x) 
calculation was performed in The general method ^^'^^ of structure function 
calculations in the QCD sum rule approach was used. In this method the structure 
function at intermediate x can be found. The main difference in comparison with 
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is that in the case of hi (x) determination the chirahty violating structure is studied. 
As a result hi{x) was found to be proportional to the quark condensate < | mm | > 
with the correction term in OPE proportional to the mixed quark-gluon condensate 
g < \ uan^G^^^,\"'u | >. It was obtained that for proton h^^x) S> hf{x), and as a 
consequence hi{x) ^ (4/9)/?"(,t). The calculations in are valid at 0.3 < x < 0.55. 
The extrapolation in the region of small x can be performed using Regge behavior 
hi{x) ~ the extrapolation in the region of large x, using the inequality (98). 

The final result of calculation of hi{x) with the extrapolations is shown in Fig. 5. 



Fig. 5 

The ?i-quark contribution to the proton 
structure function hi{x) based on QCD sum 
rule calculation at intermediate x, 0.3 < 
x < 0.55. At a; < 0.3 an extrapolation ac- 
cording to Regge behavior was performed, at 
X > 0.55 (dashed line) the saturation of in- 
equality (98) was assumed. 



It is expected that the accuracy of /i"(a;) determination is about 30% at a; ~ 0.4 and 
about 50% at X — 0.6. The inequahty hi{x) > gi{x) suggested in was confirmed. 
Numerically hi{x) is rather large, what gives a good chance for its experimental 
study. 

It is difficult to produce intense high energy polarized antiproton beam {pp col- 
lisions are preferable for the study of hi structure function, since the same enters 
(102) for target and projectile). For this reason there are suggestions of experiments, 
where only target proton is polarized, but also the polarization of the final particle 
is measured 

12. Conclusions 

In the last year the situation, connected with the so called "proton spin crises" be- 
came less dramatic, than it was before, after the appearance of the first EMC data. 
Although there is no contradiction between the new and the old data, the world aver- 
ages moved towards theoretical expectations. The simplest theoretical prediction — 
the EUis-Jaffe sum rule value moved slightly towards experimental number after 
appearance of new data on hyperon /9-decay. As a result, the gap between the 
experimental value of Tp and EUis-Jaffe sum prediction shrinked. The new neutron 
data are more or less in agreement with theoretical expectations and there is no 
contradiction with the Bjorken sum rule. 
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On the other side the theoretical understanding of the nucleon spin content 
increased essentially: we understand now that EUis-Jaffe sum rule is not the last 
word in the problem, we have more knowledge about the role of gluons and strange 
quarks. Since the gap between the theoretical expectation given by Ellis-Jaffe sum 
rule and world averaged experimental value became narrower, a possibility appeared 
to share this difference between the sourses, nonaccounted in EUis-Jaffe sum rule 
(strange sea, gluons, higher twist terms) in more or less conservative way. For 
example, such a scenario is not in contradiction with experimental data: at — 
10 GeV'^ Ag ^ 0.5 - 1.0, As ^ -0.02 - 0.05, E f« 0.4 - 0.5 and higher twist 
corrections comprising ~ 4% (the latter two times less than in Table 1). 

Evidently, these numbers would surprise nobody and would not require any 
revolutionary changes in our understanding of nucleon structure. The problem is: if 
this or a similar scenario is realized in the nature. Future, more precised experiments 
on deep inelastic /x(e) nucleon scattering will help us to obtain the answer this 
question. I emphasize, that these experiments, even very precised, are not enough 
for the final solution of the problem: an independent experiments, measuring the 
part of the nucleon spin carried by strange quarks and gluons are needed. A new era 
in the investigation of the nucleun spin structure will open, when the measurements 
with the transversally (relative to the beam) polarized nucleon will start. May be 
new surpises are waiting us here. 

I am grateful to J.Lichtenstadt for the discussion of experimental results and 
for the information about new data and to M.Karliner for helpful discussion of new 
theoretical developments. 
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